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Abstract 



We consider non-renormalizable 1/Mx interaction terms as a perturbation of the conventional 
neutrino mass matrix. Particular attention is given to the gravitational interaction with Mj = 
Mpi. We find that for the degenerate neutrino mass spectrum, the considered perturbation 
generates a non-zero U e 3 which is within reach of the high performance neutrino factories and 
just on the borderline to be of interest for supernova physics. For the hierarchical mass spectrum 
this effect is small. For 1/Mx interaction terms with Mx about the GUT scale, a detectable 
Z7 e 3 term is induced for the hierarchical mass spectra also. Numerical estimates are given for all 
the above mentioned cases and renormalization effects are considered. 

Introduction 

One of most important issues in neutrino physics is the magnitude of U e 3 or equivalently of the 
mixing angle #13. The actual value of this parameter is of great interest to various aspects of 
neutrino physics: To theory, to the search for CP violation in the next generation of long baseline 
oscillation experiments, to interpret supernova neutrino signals, etc. The only solid information 
we have on this parameter is the upper bound that is based on the CHOOZ experiment [Q, that 
is $13 < 7° at 1 sigma (13° at 3 sigma). 

In this paper, we point out that a tentative lower bound on this parameter arises from physics 
above the grand unified theory (GUT) scale, in particular from the gravitational interaction of 
neutrinos. Let us describe this idea. Most probably, grand unified dynamics generates the 
main part of the neutrino mass matrix. However, contributions from other sources are likely to 
exist. Specifically, gravitational interactions [21 El H] can produce additional terms in the mass 
of neutrinos, and these can affect the size of #13. 

The relevant gravitational dimension-5 operator for the spinor SU(2)l isodoublets, 1 ip a = 
{u a ,£ a ) and the scalar one, ip = (ip + ,p ), can be written with the operators introduced by 
Weinberg [3] as 



where M P1 = 1.2 x 10 19 GeV is the Planck mass, and A is a number 0(1). In eq.ljlj). all 
indices are explicitly shown: the Lorentz indices a,b = 1,2,3,4 are contracted with the charge 
conjugation matrix C, the SU(2)l isospin indices A, B, C, D = 1, 2 are contracted with e = ia2', 
a m (to = 1,2,3) are the Pauli matrices. After spontaneous electroweak symmetry breaking, the 

1 Here and everywhere below we use Greek letters a, /?,... for the flavor states and Latin letters i,j,k... for the 
mass states. 
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Lagrangian Q generates additional terms of neutrino mass: C mass = Xv 2 /Mp\ u a C~ x up, where 
v=17A GeV is the vev of electroweak symmetry breaking. We assume that the gravitational 
interaction is "flavour blind", i.e. A does not contain a, (3 indices. In tis case, the contribution 
to the neutrino mass matrix is of the order of: 




where the scale [i ('scale of perturbation') is 

fi = v 2 /M Pl = 2.5 x 1(T 6 eV. (3) 

In our calculations, we take eq.@ as a 'contribution of perturbation' to the main part of neutrino 
mass matrix, that is generated by GUT. In other words, our results are normalized to the case 
A = 1/2; more discussion on A follows. 

However, there might be some other interactions of the form Q with a scale Mx less than 
Mpi, e.g. from string compactification, non-perturbative dynamics, flavor physics or simply 
higher GUT scales. These interactions are not necessarily flavor blind. This implies not only 
a different scale of perturbation n = v 2 /Mx, but also a modified structure for the matrix of 
perturbation, that we denote by X a p henceforth. In the case of quantum gravity, \ a p = A. 

In this paper, we evaluate the effects of such operators on the conventional neutrino mass 
matrix, given for example by the seesaw mechanism [5J H3 IE]- We demonstrate that the 
gravitational perturbation generates a non-zero U e % even if the unperturbed mass matrix has 
U e 3 = 0. This non-zero value of C/ e 3 can be considered as a lower bound imposed by gravitational 
effects and can be translated into a lower bound for the angle #13. The operators with the 
scale Mx < Mp\, but higher than the GUT scale (string scales etc.) provides much stronger 
effects which are of potential interest to present and future observations. For reference and 
for definiteness we recall that the supersymmetric unification scale is Mqut ~ 2 x 10 16 GeV. 
However the considerations outlined below are not tied to any particular value of Mqut- 

Before passing to the calculations, we discuss in more detail the coefficient A in eq.(|T|). 
The operator in eq.Q is a term of the effective Lagrangian produced in quantum gravity, 
thus the coefficient A could be calculated if the details of this theory were fixed; but this is 
not the case at present. There is another way to see that the coefficient A is defined with 
considerable uncertainty. Indeed, an SU{2) Fiertz transformation yields (t/t a2(fip)((p t a^d <~p) = 
—2(ip <r 2 </?) (ijj &2 </?)• Thus, eq.(0) can be recast in the form 

£ g rav = ~ x ^p^" 6 ^) C '~ 1 (V^e^V 3 ) + h.c, (4) 

(suppressing all indices except those in flavor space). Here, the reader readily recognizes the 
operator introduced in Ref.[3], but with a factor —1/2 in front of it. In other words, if we 
knew that quantum gravity yielded the first operator with a coefficient equal to 1, the second 
one would be generated with a coefficient —1/2, not 1. In summary, we have to live with an 
uncertainty in the coefficient A, and this in turn reflects into an uncertainty 0(1) in our results 
for U e3 . 
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Perturbative expansion 



A natural assumption is that the unperturbed (O'^-order) mass matrix A4, 

( U el U e2 Ues \ 

M = U* diag(M i ) U\ where U= \ U„i U„ 2 (5) 




(where U a i is the usual mixing matrix, and M, the neutrino masses) is generated by grand 
unified dynamics. 2 Most of the parameters related to neutrino oscillations are known, the 
major exception is given by the mixing element U e z- We adopt the usual parameterization: 
\U e 2/U e i\ = tan #12, | £^2/^3 1 = tan #23 and |C7 e3 | = sin #13, or equivalently #12 = lj, #23 = i> 
and #13 = 4>. Note that in our approach Mi are real and non-negative and we include all possible 
phases in the mixing matrix: 

U = diag(e^) R(9 23 ) A R(9 13 ) A* R(9 12 ) diag(e^)- (6) 

The phase 5 appearing in A = diag(e^/ 2 ,l,e- i ' 5 / 2 ) is the one that affects oscillations, are 
the so called Majorana phases and fi are usually considered as a part of the definition of the 
neutrino field. It is possible to rotate away the phases fi, if the mass matrix (jSJ) is the complete 
mass matrix. However, since we are going to add another contribution to this mass matrix, the 
phases fi of the zeroth order mass matrix have an impact on the complete mass matrix and 
thus must be retained. By the same token, the Majorana phases which are usually redundant 
for oscillations have a dynamical role to play now. 

Non-GUT effects related to a larger mass scale Mx > Mqut will add other contributions 
to the mass matrix and in particular will affect the magnitude of U e3 . Thus, let us assume that 
the mass matrix is modified as: 

M -> M + fi A, (7) 

2 

with {i = -jj^— and A being a matrix of dimensionless terms as discussed in the introduction 
PJ). The impact of the new terms on the mixing can be seen by forming the hermitian matrix 
(A4 + (j,X)<(A4 + //A), which is the matrix relevant for oscillation physics. To first order in the 
small parameter fi, the above matrix is «M' Ai + fiX^M + Ai^fiX. Now by using eq.© and the 
fact that this new mass squared matrix must be diagonalized by a new mixing matrix resulting 
in corrected eigenvalues, one can write, 

U(M 2 + m)M + Mm)U ] = U'M ,2 U'\ with m = n U f A U. (8) 

Here M and M' are the diagonal matrices with neutrino masses at th and at 1 st order in /i. It 
is clear from eq.® that the new mixing matrix can be written as 

U' = U (1+iSd), (9) 



For instance, the complete seesaw formula in minimal SO (10) gO is: M = jr (Y u Y~ x Y*+$, Y) where M x = V 
is the vev of the singlet contained in the 126-plet, of the order of or somewhat smaller than the supersymmetric 
grand unification scale. Here Y v is the neutrino Yukawa coupling, Y the 126-plet coupling to the fermions and 
£ v /V is the effective vev of the triplet in the 126-plet (these two give the noncanonical, or type II seesaw). 
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where 59 is a hermitian matrix that appears at first order in \x. 
From eq.fJHJ one obtains 

M 2 + m)M + Mm = M' 2 + [iS9, M >2 ] . (10) 
Therefore to first order in /i, the mass squared differences AM| = M 2 — M 2 get modified as: 

AM- 2 = AM 2 j + 2 {MiRe[mn} - M ,\W in,/) (11) 
and the new contributions to the mixing matrix are: 



iReimij^Mi + Mj) hn[ mij }(Mi - M j} 
3 AM 2 AM 2 



The diagonal elements of 59 a are undetermined, as follows from the phase invariance of eq.(|5]l. 
Thus we set them to zero. Putting together ea. (|12j) and the definition of m in eq.((5j). we obtain 
the contribution to U e3 at 0(h): 

5U e3 = fi(M 3 + M^/AM'i x U el Re([/'Al0i 3 (13) 
-i n(M 3 - Mi)/AM' 3 2 x U e i Im(i7 t AC/)i3 + (1 -» 2). 

One observes that it is the atmospheric neutrino mass difference AM%\ that enters into the 
final expression for 5U e3 . The solar neutrino mass difference has no role in determining the 
magnitude of this effect. If the corrected mass squared difference AM'%\ is almost the same as 
the original mass squared difference AM 2 -^ (as is true for Planck scale effects), then the above 
formula simplifies to the one given below: 



6U * = » U ^ [ m 3 _ Mi ~ l M 3 + Ml ) 



+ (1^2). (14) 



Ea. (|13j) should be used to calculate the contribution for scales less than the Planck scale. 
Eas. (|13|) and (|14|) are our main results. They describe the contribution to U e3 coming from 
a perturbation above the GUT scale. 

From eas. (|ll|) and (|12[) one can also obtain the size of the deviations of #12 and #23 from maximal 
values and study the stability of the AM 2 - under effects due to physics above the grand unified 
scale. In this paper however, we focus on U e3 . 



Numerical results 

To estimate 5U e3 numerically from ea. (|14jl we need to know the mixing terms U a j, the mass 
squared differences AM| and the absolute neutrino masses Mi, M2 and M 3 . While the former 
two can be taken either from experimental data (except for the phases that at present remain un- 
known) or from a specified theoretical model, the 'absolute' neutrino masses cannot be obtained 
from oscillation experiments. See for instance We take the solar neutrino mass difference 
AMfj = 7.1 x 10~ 5 eV 2 and the solar mixing angle 612 = 34°. The atmospheric neutrino mass 
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difference is taken to be AM| X = 2.8 x 10 3 eV 2 and the mixing angle #23 to be 45°. The usual 
CP violating phase can be taken as 5 = 0. 

To get a numerical estimate for 5U e ^ we should consider what we know on 'absolute' neutrino 
masses. More precisely there £1X6 £1 certain number of masses: 

1) The mass of the lightest neutrino M, equal to either Mi or to M3 for the normal and the 
inverted hierarchy respectively. 

2) The "cosmological" mass M cosm = M1 + M2 + M3, which can be determined from the distribu- 
tion of matter on large-scales and the anisotropy of the cosmic microwave background radiation. 
From the first paper in [Hj, we quote M cosm < 0.71 eV at 95 % CL. This bound is valid under 
the simplest cosmological assumptions and depends on them rather crucially. See e.g. [TUj . 

3) The effective mass of the electron neutrino coming from the tritium experiment M 2 = 
£ < (2.2 eV) 2 at 95 % CL [II]. 

4) Neutrinoless double beta decay constrains m ee = \ J2i U^M^ < 0.38 h eV at 95 % CL ^2] 
(h = 0.6 — 2.8 quantifies the uncertainty in nuclear matrix elements). If the Majorana phases 
are the same, e.g. zero, there is no cancellation of the 3 contributions and the bound implies 
Mi f« M Ue « M cosm /3 < 0.38 h eV. If the Majorana phases produce the largest possible cancel- 
lation, this relaxes to < 1.2 h eV |13j . 

The introduction of these unknown (although constrained) quantities results in uncertainties 
in the values of U e s that we obtain in our framework. In particular, we see that the second term 
between brackets in ea. (|14(l decreases when the scale M increases, while the first term becomes 
instead larger. In other words the contribution of new physics at the scale Mx to U e s can be 
rather large in the case of degenerate neutrinos, although it is possible to diminish this effect 
by certain choices of the phases. It should be noted that the relative phases between the GUT 
contribution (phases fa) and the perturbative contribution influence the magnitude of U e s. 

Planck scale effects 

When we focus on Planck scale effects, we assume that A Q/ g = A for each a and (3. Let us 
set in the unperturbed mass matrix the mixing #13 = and also fi = ai = 0. From ea.([14|) 
one observes that the surviving contribution to C/ e 3 is given by the first term. All our results 
can be easily modified to include a specific GUT model for the unperturbed mass matrix A4, or 
for different experimental inputs and phases. The dimensional factors fJ,(M% + Mi)/AM 3 2 1 and 
^(M3+M2)/AM| 2 are equal with good accuracy and they are multiplied by cos #23+sin#23 ~ v2 
|14j . Note also that the dependence on the solar mixing angle in this case disappears. Finally 
from ea. (|14j) we get the range of values for U e s : 

^3|pi anck = 7xHr 5 -6xlO- 3 , (15) 

where the lower limit corresponds to the normal hierarchy for Mi = (AM| 2 ~ AMjj = 
2.8 x 10~ 3 eV 2 \14\ ) and the upper limit corresponds to the kinematical limit from tritium 
decay search, M < 2.2 eV. This corresponds to what is usually called the "quasi degenerate" 
spectrum, where the common mass scale is much higher than the splittings between the masses. 

Some remarks are in order: 
(i) For the case of the inverted hierarchy, the numbers are practically the same as that for the 
normal hierarchy, {it) If one takes the cosmological limit on neutrino mass, M < 0.71 eV into 
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Figure 1: Impact of the Majorana phases on calculated sin 2 2#i3 values for Planck scale effects 
and a degenerate (tritium mass) spectrum, for ai,a 2 = — 180°. For illustration we show the 
contours where the value is: 1 • 10~~ 5 , dotted line; 4- 10 -5 , dashed line; 1 • 10 -4 , long dashed line. 
The value at the peak (innermost region) is 1 x 10~ 4 . The value at the lower left corner (where 
ai = a 2 = 0) is 3 x 10~ 6 . 



account, then the upper limit becomes 5 x 10 4 , which is a order of magnitude lower than the 
one in ea. Q15JI . (Hi) The range in ea, ()15|) can be alternatively written as: 

sin 2 20i3 = 2 x 10~ 8 - 1 x 10" 4 , or as 13 = 0.004° - 0.3°. (16) 

(iv) When the cosmological constraint is used, the upper limit of sin 2 2#i3 is 1 x 10~ 6 , corre- 
sponding to an angle #13 of 0.03°. 

These numbers can be regarded as lower limits to U e ?, from Planck scale physics. They are 
admittedly small for a hierarchical spectrum, but large for degenerate neutrinos. 

The impact of phases 

For arbitrary values of the left and the right phases these numbers will change. Note that 
<23 can always be taken as 0. Putting arbitrary values of phases for the hierarchical spectrum 
gives values ranging between 1.6 — 1.9 x 10~ 8 for sin 2 2^13. So one observes that the phases have 
only a mild influence on the value of U e %. 

For a degenerate spectrum, 3 phases can affect the bound significantly. For illustration in 
fig.© we give sin 2 26\^ as a function of the right phases a\ and a 2 as a contour plot, for fixed 
arbitrary values of /1 = 66°, f 2 = 20° and f^ = 48°. Observe that the maximum value of 
1 x 10 -4 given in ea. (J16|) is reached for some combination of the phases. We also see that some 
particular choice of phases can heavily suppress sin 2 2#i3 and make it vanishingly small. 

3 For degenerate neutrinos the value of U e 3 usually scales as/jx M. So the results in fig.Q can be easily scaled 
to other values of M. This applies also when fi > v /Mpi — see next item. 
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Table 1: Impact on U e3 of a flavor blind mass matrix of perturbation at the scale Mx = 10 
GeV (left) and at the scale Mx = 10 17 GeV (right), for the two types of mass spectra (first 
column). 

Effect of scales below Mp\ 

The lower bounds we obtained increase if the scale of new physics decreases (as it is clear from 
ea.([14|)). In other words any model that needs a scale Mx above Mqut but below the Planck 
mass, can lead to a potentially interesting contribution to the still unknown mixing parameter 
U e3 . In this manner the scale Mx can lead to a significant U e3 even for non- degenerate spectra, 
in contrast to the case for Planck scale effects. 

We illustrate this in table (JTJ for two scales, each successively one order of magnitude less 
than the Planck scale. Here we set again all the phases to zero. We would also like to point 
out that for scales lower than the Planck scale and for a degenerate neutrino spectrum, the 
corrections to the solar neutrino mass difference become large (as can be seen from ea.(|ll|>) and 
fine tuning of the phases becomes unavoidable to control the corrections. This can be taken to 
suggest that for scales below the Planck scale, the degenerate spectrum is unnatural for such 
effects. Equivalently one can also say that if the degenerate spectrum is hinted at by some other 
phenomena, then it strongly constrains the scale of such flavor blind contributions. 



Comparison with renormalization effects 

The neutrino mass matrix is renormalized from the GUT scale to the electroweak scale Mz ~ 
Mhiggs, already due to the couplings of the standard model. This renormalization modifies 
U e3 and possibly contributes to make it non-zero. The purpose of this section is to evaluate 
quantitatively this effect. 

At one-loop, renormalization in the standard model is described by, 



M Q/ 3 -> rj a M a /3 V(3, with rj a = fj exp 



MauT ?(n\^9. 

(47TJ Z JMz Q 



(17) 



The indices a, f3 should be not summed over. The coefficient 77 (mostly due to gauge and top 
Yukawa couplings) modifies the 3 masses, i.e. Mi — > r/ 2 Mj, but leaves the texture of the mixing 
matrix untouched since it is flavor blind. The other contribution is different for the 3 flavors. 
Thus it changes the texture and the mixing angles. It arises since the Yukawa couplings y a of 
the charged leptons are different and the numerical coefficient b is not zero. For instance, in 
the standard model b = 1/2, while in its supersymmetric version, b = —1 (here and below, 
we consider supersymmetry at the electroweak scale). The largest part of this correction is due 
to r Yukawa coupling y T . When this renormalization is small we can apply the perturbative 
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Table 2: Impact of renormalization on U e3 for the two types of mass spectra (first column) in 
the super symmetric standard model. The value of tan j3 is equal to V3 in the left part of the 
table and to 10 in the right part. 

expansion developed above. In fact, including the effect of r\ in Ai (by redefining rj 2 Ai — > A4, 
so that i] 2 Mi — > Mi and U — ► U) we get an expression similar to eq.(|7)>: 

/ 7W eT \ 
M -> M + e • AV , when e = r] T /r) - 1 « 1. (18) 

TT I 

In the basis where the unperturbed mass matrix (the first one in ea. Q18Jl ) is diagonal, the 
components of the term of perturbation (the second term) are, 

mij = e ■ [Mi U*i U rj + U ri U* T] Mj]. (19) 

The results of eqs. ([8H4j) still apply, when we replace the matrix m introduced in eq.(jHJ) with the 
one given here (this is the reason why we use the same symbol). 

Let us focus now on the case of interest when the mixing U e 3 is at the GUT scale. Hence 
we have U T \ = sin #12 sin #23e Ml , U T 2 = — cos 612 sin^e" 12 and U T 3 = cos #23 (from ea. (|17JI we 
see that the phases fa do not play any role in these considerations and can be set to zero). The 
contribution to U e3 from renormalization is: 

5U e3 = £ - ■ sin 20i2 sin 2023 (/1 - h) {^~^) 2 . ( 20 ) 

where fj = e ia > [cos aj(M 3 + Mj) 2 - i sin aj(M 3 - Mj) 2 }/ AM' 3 2 . A crucial point to be noted 
is that we use this equation taking the values of O12 and #23 from low energy data. This is 
correct when the considered renormalization is a perturbation. The limiting expression given 
above is obtained with two other simplifying assumptions: Firstly we assume a degenerate mass 
spectrum Mi ~ M2 ~ M at the GUT scale and secondly the Majorana phases a\ and 02 are set 
to 0. 

Note in passing that the contribution to the solar splitting AMfi from y T renormalization is 
2e sin 2 #23 [(M 2 + Mf ) cos 26i 2 + AMfi] . Thus the smallest splitting will get a contribution from 
the absolute neutrino mass M ^B] j unless #12 = 45° ^7] . This value of #12 is however disfavored 
at ~ 3<r by combined analyses of solar and reactor neutrino data ^HI- This suggests a 'natural- 
ness' criterion, viz the correction should not exceed AM|i, that means: e < 0(1) (0.01 eV/M) 2 . 

In order to describe the impact of renormalization on U e3 , the key quantity is e, given in 
ea.l|18|). To be specific, we consider the case of the super symmetric model where the grand 
unification program is most commonly implemented for a variety of reasons. The tau Yukawa 
coupling is given by y r = m T / (v cos 0) and we take the range f3 = 60° — 89°, to have perturbative 
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Yukawa couplings till Mqut (note that this Yukawa coupling is always larger than the standard 
model one, m T /v ~ 10~ 2 ). A numerical integration of the one-loop system of equations as in 
QH yields: 

- e = 7 • 10" 5 - 0.15, when (3 = 60° - 89°. (21) 

The naive estimate -e ~ y 2 T (M z ) ■ \og(M GUT / 'M z ) / '(4vr) 2 would give 1CT 4 - 7 • 1CT 2 . 

At this point we have all the ingredients to evaluate the contribution to U e $ from ea.(|2U|). 
The result is given in table (J2J) for two values of /?. This shows that renormalization effects are 
not necessarily large, even though they lead us to expect that U e $ is non-zero at the electroweak 
scale. Some remarks are in order: 

(a) The Majorana phases affect the result. For simplicity and also for the sake of argument, we 
focused the discussion on the case a\ = 02 = 0. 

(b) The result depends rather dramatically on the unknown parameter /?. 

(c) The naturalness criterion on the absolute mass scale mentioned above is satisfied for the 
lowest values of (3 and also for the largest value of the common neutrino mass considered. 

(d) The largest values of/? and of M taken together, produce inaccurate estimates in perturbation 
theory and should be considered only for illustration. 

Discussion 

A priori, there is no strong reason to believe that #13 is exactly zero. In the literature there are 
a number of theoretical models for U e %. Some of them have large values, close to the present 
CHOOZ bound (for instance some minimal SO (10) models |22U23| or models with U(l) selection 
rules [2IJ where #13 ~ 0q - @c = 13° being the Cabibbo angle). More possibilities are 
reviewed in [20] , and there are also works where U e s is generated at the weak scale via radiative 
corrections to some texture defined at the high scale |21j . 

However, the experimental indications of an almost maximal mixing angle (i.e. #23 ~ 45°), 
could suggest the view that the mixing angles take very special values and perhaps #13 is really 
very small. More generally we believe that it is appealing to speculate on the possibility that 
U e 3 has an anomalously small value, say $13 < 9 C ~ 3°. In the present paper we have argued 
that such a small value might be a window for physics above the grand unification scale. In this 
sense even a negative result from future search of #13 could be of great interest. 

Our analysis concerns the aforementioned case. We assumed that U e $ is zero at the GUT scale 
and we addressed the question on whether the physics above the GUT scale can be responsible 
for a non-zero value of The outcome is not discouraging, especially if neutrinos are mass 
degenerate, or if the neutrino masses receive contributions from other scales below Mpi- In fact, 
a sizable part of the range of values of #13 obtained here is within reach of the high performance 
neutrino factories. See for instance 2fi\. The largest values of #13 can have an effect on supernova 
neutrino fluxes via the MSW effect |27j . Indeed, the 'flip probability' in the supernova mantle 
due to the atmospheric neutrino AM 2 is approximatively given by Pf « exp[— £ (U£ 3 /10 )] 
where £ = (15 MeV / E u ) 2 / 3 . The tentative lower limit we discussed suggest that #13 could well 
be at the border of the region of adiabaticity and hence can lead to a distortion in the spectrum 
of supernova neutrinos 28 . 
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